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A NONLINEAR THEORY FOR SHELLS WITH SLOWLY 
VARYING THICKNESS 

MARTA LEWICKA, MARIA GIOVANNA MORA AND MOHAMMAD REZA PAKZAD 



Abstract. We study the r-liinit of 3d nonlinear elasticity for shells of small, 
in ' variable thickness, around an arbitrary smooth 2d surface. 

Oh: 



1. Introduction 



The following question receives large attention in the current literature on elas- 
r^ . ticity [T]: which theories of thin objects (rods, plates, shells) are predicted by the 

Mh' 3d nonlinear theory? For plates, this problem has been extensively studied through 

r~| . F- convergence; first by LeDret and Raoult [5] , leading to a rigorous derivation of the 

membrane theory, and later by Friesecke, James and Miiller [4], for the Kirchhoff, 
von Karman and linear theories. In this framework, much less has been done for 
shells. The membrane and the bending theories were obtained in [5] and [5], respec- 
tively. More recently, the generalized von Karman and linear theories have been 
rigorously introduced and justified by the authors in |7j. In this paper, we present 
these last new results, in an extended version for shells with variable thickness and 
mid-surface. 



2. The setting of the problem and the related elastic energy 

functionals 



(N 

^' 

^^ I The three-dimensional problem. Let 5' be a 2d compact, connected, oriented 

QP ■ surface in M'^, of class C^'^, whose boundary is the union of finitely many (possibly 

none) lipschitz curves. Let n{x) be the unit normal vector, T^S the tangent space, 
and n(a;) = Vn(a;) the shape operator on S, at a given x ^ S. 

For given positive lipschitz functions gi, g2 ■ S — > (0, oo), small in L°° , consider 
a family of thin shells: 

S'^ = {z = X + tn{x); x e S, ~hgi{x) <t < hg2{x)}, < h < 1. 

The elastic energy (scaled per unit thickness) of a deformation u^ £ W^''^{S'^ ,M.^) 
is given by: 

E'^iu'') = \ I W{Vu''), 

where the stored-energy density W : M.^^^ — > [0,oo] is C^ in a neighborhood of 
SO{3), and satisfies: 

VF e M^^^ Y^ g 5-0,^3) |4^(^) ^ 0, W{RF) = T4^(F), 

WiF)> Cdist^ {F,SO{3)). 
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The two-dimensional problem. We now introduce the von Kdrmdn functional 
onS: 

(2.2) 

I{V,Btan) = 2 / (31+52)22 (x,Btan - '^{A'^)tan " -Syni(^V((<72 " 5l )")) 
+ 2^ I {91+92? Q2 [X, {V{An) - AI{)tan) , 

defined for V £ V and Btan G B. Tlie space V consists of infinitesimal isometries 
V € W'^''^{S,M.^), tliat is tliese vector fields V for whom there exists a matrix field 

Ae W^'^iS,M.^''^) so that: 

drVix) = A{x)t and A{xf = ~A{x) Va.e. xeS Vr e T^S. 
The finite strain space B (see [7]), consists of the following symmetric matrix fields: 

B={l^ - lim symVu;''; w^ e M/^^'^(5,M3)). 

In (j2.2p K > is a parameter, and the positive definite quadratic forms 22(2;, •) are 
defined as follows: 

Q2{X, Ftan) = min{Q3(i^); (F - F)tan = 0}, QsiF) =. D^Wild){F, F). 

The quadratic form Q3 is defined for all F e K'^^^, while Q2{x,-), for a given 
X G S is defined on tangential minors Ftan of such matrices. Recall that the 
tangent space to 5*0(3) at Id is so(3). As a consequence, both forms depend only 
on the symmetric parts of their arguments and are positive definite on the space of 
symmetric matrices [3 . 

Relative stretching and bending. The two terms in ()2.2p are strictly tied to 
the deformations of the geometric mid-surface of 5''. Given V G V and a field 
w G W^''^{S,M.^), consider the deformations: 

4,'' = id + 1/2/1(32 - 9i)n, (/)'' ^4>'' + hV + h^w. 

Then (j)^{S) is the geometric mid-surface. A straightforward calculation shows that 
for any r e T^S: 

|a^0''|2_|a^0''|2 =2/iV^(symVu;- 1/2^2 -l/2sym(AV((52-5i)n)))r + C'(/i^). 

Hence (putting n ~ \ and Btan = symVw), the expression in the argument of Q2 
in the first term of (12. 2|) describes stretching, namely the second order in h change 
of the first fundamental form of 4)^{S). 

The second term of (J2.2I) relates to bending, that is the first order in h change in 
the second fundamental form of 4''''{S). To see this, let II'' be the shape operator 
on (t)'^{S). Indeed, for t G T^S one obtains: 

(V(/)'')-in''(V0'')T - {V4>'')-^U{V(t>'')T = h(driAn) - AUr) + Oih^). 

Notice also, that the term — l/2sym(AV((52 — 5i)"')) which is new with respect 
to analysis in [7], disappears when gi ~ 32 or cquivalently S ~ 0''(5). This term 
expresses the first order in h deficit for V from being an infinitesimal isometry on 
(j)^{S). It is because for any 77 = dr(p{x) where r G T^S we have: 

d^{V o (0")-i) ■r^ = drV- dr^'' = -/i/2T^sym(AV((.92 - gi)n))T. 
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3. The r-CONVERGENCE OF THE ELASTIC ENERGY FUNCTIONALS 

Theorem 3.1. Let e be a sequence of positive numbers, for which we assume: 

(3.1) lim e'^/h'^ = k^ < oo. 

(a) For any sequence of deformations u^ € iy^^^(S''', R'^) satisfying: 

(3.2) E'\u'') < Ce\ 

there exist a sequence Q^ £ SO{i) and c'' S M'^ such that for the normalized rescaled 
deformations: y^{x + tfi) = {Q^Y' u^ {x + htn) ~ c^ defined on the common domain 
S^ , the following holds. 

(i) y'' converge in H^^'^(5^) to the projection n : S^ > S along the normal 

n. 
(ii) The related scaled average displacements: 

{V%'']){x) = —= -f y'\x + tn) -{x + htn) dt 

Ve J-gi(x) 

converge (up to a subsequence) in W^''^{S) to some V" £ V. 
(iii) j^syuiW'^ly'^] converge (up to a subsequence) weakly in L'^{S) to some 
Btan e B. 

(iv) liminf;,^o JKE'^iu'^) > I{V,Btan)- 

(b) Conversely, for every V ^ V and Btan G B, there exists a sequence u^ E 
W^'^(S''',R^) satisfying i3.S\) and such that (i), (ii), (iii) hold (with Q^ = Id and 
c'' ~ 0) and moreover: \m\h^o -^E^{y^) = I{V,Btan)- 



Proof. Here we outline the main steps of the proof of Theorem 13.11 We refer the 
readers to [7] for a detailed exposition in the case gi = g-i = const. All convergences 
below are up to a subsequence. 

1. A first major ingredient is approximating u^'- by W^'^ matrix fields R*^ : S -^ 
SO{3), thanks to the scaling invariant nonlinear quantitative rigidity estimate from 
[3j. Noting (|3.2p and the uniform equivalence of the functionals Jdist^(V-, 50(3)) 
and £'''(•) on the uniformly lipschitz domains 5'', it follows that: 

\\W - R^Mlms-) < Ch^/^v^, \\WR''\\l2(s) + \\R" - Q'^Wms) < Ch-^V^ 

for some Q'' e 50(3). Further, for some skew-symmetric matrix field A E W^''^{S, M"^) 
one has: 

h/V^{{Q''fR'' - Id) — > A weakly in W^'^{S), 
h^/e''sym{{Q'YR^ - Id) — > 1/2A'^ strongly in L^{S). 

2. Define now: 

w''{x + tfi) = h/y/^{iQ''fu'\x + htri) - c^ - {x + htn)), 

with c'* such that Jgi w^ — 0. To prove (i) and (ii), one first uses the above 
convergences to obtain that: 

Vtanw'* — > An{Id + tn)-^ and dnw^' — > 0, strongly in L^{S^). 
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By the Poincare inequality, w^ must now converge to V^ o vr, strongly in W^'^{S^). 
For (iii), one equates various terms in ^symVy''[i/''] to notice that it converges 
weakly in L'^{S) to: 

(3.3) j {symG)tandt+'^{A'^)tan + ^sym{A\7{{g2-gi)n))=:Btan. 

where G is the weak L'^{S^) limit of the matrix fields: 

G'*(x + tn) = -^{{R''fWu''{x + htn) ~ Id). 



The convergence of G'' [7J, which is the ve'* order term in the expansion of the 
nonlinear strain (Vu'')"^Vm'' at Id, plays a major role for the expansion of the 
energy density: W{Vu'') = Vl^(Id + Ve^G''). 

3. Towards proving (iv), the previously established convergences and the definition 
of Q2 imply: 

(3.4) 

liminf -^E'\u'') > - / / Q2{x, (sym G)tan) dtdx 

22(2;, (sym Ga)tan + (s + (52 - gi)l2)(VtanA)ri\ dsdx, 



1 



(91+92)72 



2 JSJ-(gi+g2)/2 

where we used the fact that: 

(sym G)tan{x + tn) = (sym Go)tan(a;) + t{VtanA)n 

is linear in t. Calculating (sym G())tan from (j3.3[) and collecting various terms in 
the argument of the quadratic form Q2 above, we see that the right hand side in 
()3.4p equals the von Karman functional /(V, Btan)- 

4. If V^ e V n W^-'°°{S,M?) and Bta„ = symVw with w € W^'°°{S,R^), then the 
recovery sequence is given by the formula below. When V and Btan do not have 
the required regularity, one proceeds by approximation and truncation, as in [7]. 



y''{x + tn) =x + h/2{g2 - gi)n{x) + ^/e^/hV(x) + Ve^w{x) 

+ h{t - 1/2(52 - gi))n{x) + {t- 1/2(52 - gi))V?^(nVtan - V(Fn)) {x) 

- h{t - 1/2(52 - 5i))\/^n^Vw + (t - 1/2(52 - gi))hV^d%x) 

+ l/2{t - 1/2(52 - gi))^hV^d\x). 

The vector fields dP,d^ G W^'°°{S,M.^) are given by means of the linear map 
c(x,Ftan) which returns the unique vector c satisfying: 

Q2ix,Ftan) = min{Q3(Fta„ + c(g)n(a:)+n(a;)(g)c); c G M^}. 



Then: 



d° =2c (x,Btan - '^iA%an ~ -syui {AV {{g2 - gi)n))j 

+ nA^ri - '^{ffA^n)n + -(^AV{{g2 - gi)n)^ n 
d^ =2c(a;,sym {V{An) ~ AlVjtan) + n^ AI{- ffV{An). 
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4. The dead loads and convergence of minimizers 

Consider a sequence of forces /'' £ L'^{S,M.^) acting on S, with the foUowing 
properties: 

(4.1) [{91+92)1^^0 and Urn ^=/" = / weakly in 2.2(5). 

Define their extensions f'''{x + tn) = det (Id + in(a;))^ f'^ix) over S''*, and the 
maximized action of /'' over aU rotations of S'^ as: 



m ' = max 



1 



Qeso(3) h Jgh 
The total energy functional on S^ is: 



f'{z)-QzAz. 



J''{u'') = E^{u^) + Tn^ -- f f'^u^. 

h Js" 

As in [7j, one can prove that: 

> inii^j'^iy''); u^ G W^'^{S^,m?)\ > -C. 

We further introduce the relaxation function r : S0{3) — > [0, oo], with its effective 
domain TW = {Q G 50(3); r{Q) < oo}: 

(4.2) r(Q) = inf I liminf^fm''-! / f'^-Q^zdz); Q'* G 50(3), Q'^^q]. 

1 h^o e" V h J gh / J 

Theorem 4.1. Assume \3.1]) and J^.J] ). Let u^ e W^''^{S,M.^) he any minimizing 
sequence of -^J^, that is: 

(4.3) limflj''(y'')-inflj'')=0, 
h—>o \e e"- J 

Then the conclusions of Theorem \3.1\ (a) hold, and any accumulation point Q of 
{Q^} belongs toA4. Moreover, any limit {V, Btan,Q) minimizes the following func- 
tional, over V e V, Btan G B and Q £ M: 

J[V, Btan,Q) = I{V, Btan) ^ / (^1 + 92) f ' QV + r[Q). 

The proof follows as in [7] . An equivalent formulation of Theorem 14.11 in terms 
of F-convergence is possible. 

Example. When f'^ — hye^ f and gi = g2, then: 



M^ \q£ S0(3); [ f ■Qx= max f f ■ Qx 

[ Js QGSO(3) is 

(in the general case the inclusion C is still true, but not the other one). Further, 
r = on Al, so the term r{Q) may be dropped in the definition of J. In the general 
case, both r and M depend on the asymptotic behavior of the maximizers of the 
linear functions 50(3) B Q i-^ J^^ f^{z) ■ Qz dz. 

Approximately robust surfaces and higher scalings. Some classes of sur- 
faces, including surfaces of revolution, developable surfaces with no flat part, and 
elliptic surfaces, have the ability to recompense the second order stretching (of 
^^(5)) introduced by the infinitesimal isometry V, through a suitable second order 
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displacement. As an example, surfaces with flat parts do not enjoy this property 
(for any gi, §2). For such surfaces [7], described by condition: 

{A^ + sym(AV((52 - gi)n)); F G V} C 6, 

the limit elastic energy I{V,Btan) simplifies and should be replaced (in Theorems 
I3.1l and l4.1l regardless of k > 0) by: 

(4.4) i{V) = l-^j[g,+g^fQ^[x, {V{An) - AU) ^^ 

When K = and gi = (72, the stretching term is negligible for any S and hence 
the von Karman theory again reduces to the linear theory, with the elastic energy 
given in (|4.4p . For more discussion, see [7] . 
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